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Abstract. We report on the results of a two epoch study of the low power Compact Symmetric Object 4C 31.04. Observations 
performed with the VLB A at 5 GHz in 1995 and 2000 have yielded images of this source at milliarcsecond angular resolution. 
A central core is detected, with bright compact hot spots and extended lobes on both sides. Model-fitting and other analysis 
of the data (brightness profile, difference map) clearly indicate that the source is expanding. We estimate the velocity of this 
expansion to be (0.085 ± 0.016) mas/yr, i.e. (0.33 ± 0.06) c in both hot spots. Assuming a constant expansion velocity, we 
estimate the kinematic age of the source at 550 yrs. We also study the spectral index using VLB A observations at 1.3 GHz and 
MERLIN at 22 GHz. The derived spectral age is 3000-5000 years in equipartition conditions. The two estimates are discussed 
and found to be in agreement, given present uncertainties. 
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1. Introduction 

In recent years, considerable efforts have been devoted to the 
study of Compact Symmetric Objects (CSOs). The name for 
this class of objects is derived from their morphology - they 
have sub-galactic dimensions but the same symmetric, two- 
sided aspect as classical large scale FRI and FRII radio galax- 
ies. The CSOs also show other interesting properties in the 
radio: (1) their spectrum can peak at frequencies as high as 
a few GHz (hence their frequen t clas sification as Gigahertz 
Peaked Spectrum - see O'Dea (1998), Polatidis & Conway 
(|2002[) and Polatidis et al. (|2002|) for recent reviews), (2) the 



an evo l utionar y sequence (Fanti et al. |1995| ; Readhead et al. 
19964 |l996b| ; Polatidis & Conway |2002| ) in which a radio- 



rate of Hi absorption is higher than in common radio galaxies 
(e.g. Van Gorkom et al. [1989^ Pihlstrom |2001| ), and (3) vari- 
ability and polarization are almost non-existent (Fassnacht & 
Taylor 2001). At other wavelengths a comprehensive study of 



the properties of CSO host galaxies is still missing, although 
HST observations have suggested that in at least some cases the 



quiet galaxy turns into the host of a large, kpc-scale radio loud 
object. It is therefore of great interest to (1) investigate evidence 
supporting and quantifying the age of CSOs and (2) relate their 
current properties to the mechanisms triggering the nuclear ac- 
tivity. Indeed, not only the triggering but also the fueling of this 
activity can be studied in these objects. Since their orientation 
provides continuum emission on both sides of the center of ac- 
tivity one can search for either Hi or free-free absorption (Peck 
et al. |2000| ). 

The radio galaxy 4C 31.04 (B2 1 16+31) was observed 
with VLBI by Wrobel & Simon ( |1986| ) at 92 cm using 7 anten- 
nas. Their image showed two components separated by about 
70 mas. Cotton et al. (1995) tentatively classified this source 
as a low redshift (z = 0.0592) CSO. Giovannini et al. (2001, 



hereafter G01) confirmed the CSO structure and identified the 



host galaxies are not completely relaxed (Perlman et al. 2001 1 



Our current understanding of CSOs suggests that they are 
small by virtue of their youth (youth scenario) and not con- 
fined by an extremely dense medium (frustration scenario). 
If they are young then CSOs represent the very first stage of 



core with a weak flat spectrum central component. 

On the basis of spectral line VLBI observations, Conway 
(1996 - hereafter |C96 , 1999) has presented evidence for a cir- 
cumnuclear Hi disk, viewed close to edge-on, whose axis is 
that of the radio jet. The presence of a circumnuclear disk 
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which reveal a disk-like dusty feature of somewhat larger scale 
but similar orientation to the one suggested by Conway. The 
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optical nucleus is also found to have cone-like features well 
aligned with the radio axis. 

After a five year gap following the first observation at 5 
GHz presented in G01, we re-observed this source with the 



VLBA in order to look for a possible expansion. After describ- 
ing the observations and data reduction in Sect. 2, we present 
our detailed comparison of the two epochs and information on 
spectral index in Sect. 3. Some implications of the results are 
discussed in Sect. 4. Conclusions are summarized in Sect. 5. 

We assume Ho = 65 km s _1 Mpc -1 and q ( ) = 0.5 through- 
out. At the distance of 4C 31.04 (z = 0.0592) 1 mas corre- 
sponds to 1.2 pc and 1 mas/yr = 3.9 c. 

2. Observations 

2.1. 5 GHz Observations 

The first epoch observations at 5 GHz were obtained on July 
22, 1995 (1995.554), with the Very Long Baseline ArrayQ plus 
a single VLA antenna (Yl) for ~5 hours. Details about these 



observations are discussed in Appendix A of 301 

In the second epoch we obtained ~10 hrs of observing 
time with the full VLBA + Yl in full polarization mode on 
July 3, 2000 (2000.504). The observing frequency was cen- 
tered at 4.971 GHz with 4 IFs and both polarizations recorded. 
Correlation was performed in Socorro. All calibration and 
fringe fitting were done within NRAO Astronomical Image 
Processing System (AIPS): J0136+4751 was used as an am- 
plitude and leakage term calibrator, while J 175 1+0939, 3C395 
and J0423-1020 were used as calibrators for the polarization 
angle. 

For the purposes of this work it was important to have both 
datasets reduced in the same way. Thus we reanalyzed first 
epoch data, in order to have the same calibration procedure 
for both observations. We then exported the two datasets to 
DIFMAP^, which we used to perform the editing, self calibra- 
tion and model-fitting. To properly register the two images, we 
used the 1995.554 image as a starting model for the self cali- 
bration of the 2000.504 data. We chose a restoring beam of 3 
mas FWHM for both epochs. 

2.2. 1.3 GHz and 22 GHz Observations 

It is of great interest to compare the dynamic age - derived from 
the study of hot spot velocities, to the radiative age - which can 
be estimated from synchrotron losses. In order to make this 
comparison, observations at multiple wavelengths are needed. 

Data at 1340 MHz were obtained with the VLBA + Yl 
on July 7, 1995 in a project aimed to study the nature of Hi 
absorption in CSOs. Spectral line results and a continuum map 
with angular resolution of 3 mas (FWHM) were reported in 
C96 Details on data calibration and reduction will be presented 
in a forthcoming paper (Conway et al., in prep.). Here we point 



1 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under cooperative agreement with the 
National Science Foundation. 

2 DIFMAP was written by Martin Shepherd at Caltech and is part 
of the VLBI Caltech Software package 



out only that the same strategy used for the 5 GHz data was 
adopted. 

We also investigated the high energy end of the spectrum by 
considering a MERLIN observation at 22 GHz. Observations 
were made on 29th January 1998 in continuum mode with a 
bandwidth of 14 MHz, using all 5 of the MERLIN telescopes 
available at this frequency. The observations lasted 12 hrs. 
Observations were made using phase-referencing, but only low 
dynamic range images were obtained using phase-referencing 
so final images were made using self-calibration. Good fits to 
the phase and amplitude data were obtained. Final restoring 
beam was 10 mas (FWHM). 

We used the 1.3 GHz data and the second epoch 5 GHz ob- 
servation to map the spectral index. For a proper comparison 
we tried to achieve a similar (m, v) coverage at the two observ- 
ing frequencies by cutting the shortest baselines in the 20 cm 
data and the longest ones in the 6 cm data. Final images were 
then produced with the same angular resolution and cell-size. 
The source position in the two images was registered using the 
isolated, unresolved core. 

The resulting (u, v) coverages at 5 and 22 GHz are quite 
different, so it was not practical to produce a spectral index 
image between these frequencies. However, comparable short 
baselines are present, so that we do not expect missing flux in 
either image. This allowed us to derive an integrated spectrum 
for different regions of the source. 

3. Results 

3.1. Morphology 

In Figure [l] we show the final image of 4C 3 1 .04 at 5 GHz for 
the July 2000 epoch, convolved with a circular restoring beam 
of 3 mas (FWHM). The overall structure of the source is the 
same as was observed in 1995 (both at 5 and 1.3 GHz), con- 
sisting of a compact core component with hot spots and lobes 
on either side. The lower resolution (10 mas FWHM) 22 GHz 
MERLIN image shown in Fig. ^ shows the same features, ex- 
cept that the lobes look much fainter and are barely detected. 

At all frequencies the total source extent is ~100 mas, 
which corresponds to ~120 pc. Total flux density in the 5 GHz 
maps is ~ 1.25 Jy, i.e. about 90% of values measured with 



single dish telescopes (Becker et al. 1991; Gregory & Condon 
I99\b. This yields a total power of 1.14 x 10 25 W/Hz at 5 



GHz, which makes 4C 31.04 one of the weakest known CSOs 



(Polatidis et al. 2002) 



Within Difmap, we model-fit the visibility data of the two 
5 GHz observations with a set of 10 elliptical Gaussian com- 
ponents. We used large, extended components in order to pick 
up the lobes' flux density as well as compact components to 
describe the core and hot spots. We fixed the parameters of the 
largest components (size > 10 mas) and let the most compact 
ones vary; results of the fits are listed in Table |l|. Uncertainties 
for each component are computed from the ratio between its 
FWHM and SNR and are all within 0. 1 mas. The most compact 
components, except for the core, are all located in the hot spots 
regions. The western hot spot is labeled Wl, while the best fit 
to the eastern end is obtained using three different components 
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Fig.l. The July 2000 VLBA image of 4C 31.04 at 5 GHz. 
Contours are drawn at -0.4, 0.4, 0.57, 0.8, 36.2 mJy/beam 
by V2 intervals with negative contours shown dashed. The 
peak flux density is 40.2 mJy/beam. The grey scale range is 
from -0.75 to 40.2 mJy/beam. Note the compact core, sym- 
metrically placed hot spots, and the unusual 'hole' in the east- 
ern lobe. 
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Fig. 2. MERLIN image at 22 GHz restored with 10 mas 
FWHM. Contours are drawn at -4, 4, 8, 16 32, 48, 64, 96, 
128 mJy/beam. The peak flux density is 138.5 mJy/beam. The 
grey scale range is from -4.3 to 138.5 mJy/beam. 



-El, E2 (also the image peak of 40.2 mJy/beam) and E3. Given 
their high brightness, position close to the edge of the lobe and 
compactness, we will consider Wl and E2 as the two primary 
hot spots of this source. It is less likely that either El or E3 - 
both quite fainter than E2 - may be a primary hot spot; it is 
however possible that they are secondary hot spots similar to 



those found in many extended FRII sources. The innermost E3 
could also be related to the presence of a knot component in an 
underlying jet. 

Indeed, even if a classical jet of plasma is not detected, hints 
of a jet-like feature are clearly present in the enhanced bright- 
ness visible all along the line connecting the core, and in the 
east hot spot which has an elongated structure on the west side 
in the full resolution image published in 301 Giovaninni et 
al. suggest that the source is oriented close to the plane of sky 
(6 > 75°), with the west lobe the one nearest to us, in agreement 



with C96 



With the increased image fidelity obtained in this second 
epoch we can also confirm the presence of a 'hole' in the east 
lobe. This large region of lesser brightness lies close to the east 
hot spot, in a direction perpendicular to the jet. While visible in 
the first observation, it looks even dimmer and more defined in 
this second epoch. In contrast, the west lobe does not present 
very clear features. Being very extended and of low brightness, 
it cannot be mapped with high fidelity by VLBI technique; 
however, we are confident that this does not affect the com- 
pact hot spot region, which is well detected on most baselines. 
Fig. U and the 1.3 GHz data ( 296) also confirm that it is indeed 
a compact, flat-spectrum component. 

At their outermost extents, both lobes have a very sharp 
edge, running along a very wide front (~ 20 mas in the eastern 
lobe but even ~ 40 mas in the western one). The sharpness 
of this border is again in agreement with this region being the 
place where the plasma is advancing into the external medium 
as the source grows. 

3.2. Component Motions and Kinematic Age 

The first hot spot velocities in CSOs by Owsianik & Conway 
(1998) required time baselines of ~ 15 years, with VLBI ob- 



servations at 5 GHz. By going to higher frequencies (15 and 
43 GHz), Taylor et al. (2000) reduced this to ~ 5 years for a 



few strong CSO s. Cur rently 10 sources have been measured 
(Polatidis et al. 2002 ). Typical velocities are 0.2/iT.Lc (i.e. 



0.3 l/i 65 c), which due to the large distance of the sources, trans- 
lates into angular motions of typically only a small fraction of 
mas/yr (~32 /ias/yr at z = 0.1), thus requiring long time inter- 
vals and high angular resolutions. We consider observations at 
5 GHz separated by 5 years. More observations in the future 
are planned to extend the time baseline and to better constrain 
the results. 

To measure the advance velocity in 4C 31.04, we compare 
the model-fit of the dataset in our two epochs. Fig. || illus- 
trates the motion of all the compact components with arrows of 
length proportional to the magnitude of their motion. It is ap- 
parent that all the components have moved away from the core. 
In particular, the two hot spots Wl and E2 show motions almost 
in the same direction as their position angles. Furthermore, the 
increase in hot spots separation becomes even more clearly vis- 
ible when E2 and Wl are directly compared. In Fig. ^| we show 
a magnified inset around the west hot spot region with con- 
tours drawn for both epochs. In this image the 2000.504 pro- 
file has been shifted by 0.3 mas westward, so that the position 
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Table 1. Gaussian model and relative proper motions for 4C 31.04. 
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NOTE - Parameters of each Gaussian component of the model brightness distribution: S , flux density; r, 9, polar coordinates of the center of 
the component relative to an arbitrary origin, with polar angle measured from north through east; a, b, major and minor axes of the FWHM 
contour; <I>, position angle of the major axis measured from north through east; Ar, relative proper motion of the component; v, relative projected 
velocity in units of c (h 65 = H /65 km s _I Mpc~'), along the given position angle (p. a.). 



of component E2 is the same in both maps (the displacement 
needed in declination is smaller than the measurement errors). 
The motion of the peak of Wl is clear, confirming the results 
of model-fitting. 

In Fig. [| we show a brightness profile made in AIPS with 
the task SLICE. The two profiles are derived from the two 
epochs along the same slice. This slice cuts the two hot spots 
at a position angle of -82.6° (measured north to east), pass- 
ing close to the core (but not across the core itself because of 
a small bend in the jet). In the global profile one can easily 
see the two lobes on the opposite sides (the east one being the 
brightest) and the core region in the middle (see Sect. about 
core variability). The hot spot regions are magnified in the two 
insets, showing that both in the eastern and the western hot 
spot the 2000.504 profile is shifted outward. This displacement, 
though small, is clear both in the rising and in the falling part 
of the profile. We remind the reader that the region between the 
core and the west hot spot is poorly imaged, so that the features 
present in the slices in this region are not completely reliable. It 
is also possible that some of the changes are due to turbulence 
present in this back-flow region. 

The increase in the hot spot separation implies an expansion 
of the whole source. Since the lobes are too extended and faint 
to measure their expansion by model-fitting, we tried to directly 
compare the two images. The difference map shown in Fig. ^| 
was done with the AIPS task COMB. The features that appear in 
this image, in particular in the eastern lobe, are clearly related 
to the expansion of the source. The two close white and black 
spots are due to the outward shift of the peak of component E3. 



Moreover, the presence of a white stripe all around the outer 
edge of the lobe, though faint, is suggestive of an expansion 
of the whole eastern lobe in this direction. Note that a similar 
white stripe, even if less evident, is present all around the edge 
of the west lobe. 

As these last results (shown by Fig. |[ ||) are drawn from 
the images and not from the model, we feel that they lend sup- 
port to the detection of motion and to the estimate of its veloc- 
ity derived from the model-fits. With just two epochs available, 
and a moderately short time baseline, it is difficult to better 
quantify the accuracy of our measurements of component mo- 
tions. We hope to address this difficulty with future observa- 
tions. 

Model-fitting of the present data (see Table |l]) provides the 
most reliable information on the velocity. The east and west 
hot spots present quite different displacements, AW1 = (0.54 
+ 0.07) mas and AE2 = (0.31 ± 0.04) mas. It is possible that 
such differences may be due to a velocity asymmetry related to 
an inhomogeneous ISM, but it is not possible to discuss it more 
deeply with only a two epoch measurement. Since the two arms 
are about the same length and the source is near to the plane of 
sky, we will assume that the average velocity over the long term 
is the same in both directions. 

We therefore try to estimate the advance speed by averag- 
ing up the results for the two hot spots Wl and E2. This yields 
an average motion of (0.42 ± 0.08) mas in 5 years. This indi- 
cates that each hot spot is advancing at (0.085 + 0.016) mas/yr, 
i.e. (0.33 + 0.06) hzl c. The most distant component (Wl) is 
located 46.62 mas from the core, so it would have taken 548 ± 
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Fig. 3. The 2000 VLBA image of 4C 31.04 at 5 GHz with 
arrows indicating the positions and motions of components 
derived from model-fitting. The model-fit results are shown 
in Table 1. Components larger than 10 mas have been held 
fixed between the two epochs for reasons discussed in the text. 
Component motions have been magnified by a factor of 5 for 
illustration purposes. Those arrows pointing due south indicate 
that the component was held fixed. 



100 years to cover this distance at such speed. Therefore, based 
on kinematics, we can estimate the source age at ~ 550 years. 

3.3. Spectral Index Distribution and Spectral Age 

From the measured energy distribution of a population of elec- 
trons it is possible - under a few basic assumptions - to de- 
termine how long these electrons have been radiating. This has 
proved to be a useful tool in estim ating ages of c ompact ra- 
dio sources (e.g. Murgia et al. 1999 Readhead et al. 1996a ). In 
order to perform such a study, one needs to know the spectral 
index distribution and the magnetic field in the source. 

The image in Fig. [7] shows the spectral index for 4C 31.04 
between 1.3 and 4.8 GHz at a resolution of 3.0 mas (FWHM). 
The core spectrum is inverted| (a = -0.6 + 0.05). This is prob- 
ably due to self-absorption, even afte r accounting for possible 
flux density variability (see Sect. 3^5). In the east lobe, the hot 
spot region shows a flat spectrum (a ~ 0.1 ± 0.01) with a con- 
stant and uniform steepening toward the core direction, up to a 
~ 1 .0 (+ 0.02) at 20 mas from the core, just before the 'gap' be- 
tween the core and lobe. The spectral index distribution is more 
complicated in the west lobe perhaps on account of a stronger 
interaction of the radio plasma wit h the ISM, as suggested also 
by the lobe morphology (see Sect. 4.2). The region of the faint 
hot spot is still flat (a ~ 0.25 + 0.01) but the spectral index 
steepening is more irregular moving toward the core and to the 
north. The steepest values in the more external northern regions 
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Fig. 4. Overlay of the hot spot region in the west lobe (Wl). 
Contours are drawn in grey (epoch 1995.554) and black (epoch 
2000.504) at (1, 2, 4, 6, 8, 10) mJy/beam. The two images have 
been aligned on the position of the east hot spot E2. 




MilliARC SEC 



Fig. 5. A slice through the hot spots of 4C 31.04. The central 
small peak results from the northern edge of the core (the hot 
spots are not perfectly symmetric about the core). The black 
line shows the profile of the 2000.504 image, while the grey 
line represents the profile of the 1995.554 image. Shown inset 
are magnified versions of the profile at the east (E) and west 
(W) hot spots. 



Spectral index is defined according to S(v) ocv a . 
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Fig. 6. The difference image computed as S 2000.504 - S 1995.554. 
Note the change in the core flux density of ~4 mJy. One can 
also see the black+white pattern of motion in the eastern hot 
spot, and a faint ring of white along the edge of the eastern 
lobe. The "mottled" pattern in the lobes is a result of the low 
fidelity of the image in the presence of extended emission. 



are as high as 1.0-1.2. In the core-hot spot direction the region 
from 10 to 30 mas shows again a flat spectrum (0.1-0.4) sug- 
gesting a local interaction between radio plasma and the dense 
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Fig. 7. The spectral index map of 4C 3 1 .04 at 3 mas resolution 
(FWHM). The grey scale range is -0.5 to 1.2. Overlaid con- 
tours from 4.8 GHz are drawn at 0.5, 1, 1.5, 3, 5, 7, 10, 20, 30, 
40 mJy/beam. 



ISM. Imaging problems in this low brightness region could also 
play a role. 

In this high resolution image, one can see the region of the 
hole close to the east hot spot. The spectral index is very steep 
here, not only in the deep hole itself - where it is as steep as 
1.7 ± 0.04 - but also in external lobe region, where we find a 
= 1.0 + 0.04. 

In a population of electrons with an initial power-law en- 
ergy distribution, synchrotron losses produce a curvature in the 
original straight power-law radio spectrum. The curvature man- 
ifests itself at a critical frequency related to the electron age. 
The shape of the spectrum and the critical frequency crucially 
depend on the evolution of the el ectron pitch angle distribu- 
tion with time (Pacholczyk 197C). We used the Jaffe-Perola 
model which assumes a redistribution of electron pitch angles 
on time scales s hort c ompared with their radiative lifetimes 
(Jaffe & Perola, 1973 ). To derive the critical frequency from 
the spectral inde x info rmation, we used the Synage program 
(Murgia & Fanti 1996 ). Since a spectral age measurement can 
give only the time since electrons were accelerated, we avoided 
the core and hot spot regions where local turbulence and elec- 
tron re-acceleration is likely to be present, and investigated the 
steep-spectrum, relaxed regions where the 'oldest' electrons 
are present. The break frequency was found to be ~7.4 GHz 
in the east lobe at ~20 mas from the core, and between 6. 1 and 
7.2 GHz in the northern region of the west lobe. 

The MERLIN image in Fig. |^ lends strong support to this 
estimate, showing that very little of the lobe emission is de- 
tected at 22 GHz. This cannot be due to a lack of short spac- 
ings since the MERLIN has baselines down to 6 km. It follows 
that the lobes are undetectable because of spectral steepening. 
Fig. |8] shows the spectrum for the two hot spots, the lobes and 
the core. The steepening between 5 and 22 GHz in the lobes 
supports the presence of a spectral break around 5-7 GHz. 

The local magnetic field, which is also needed to derive 
the lifetime of radiating electrons, was calculated assuming 
equip artition. We note that a recent X-ray study of a relatively 
compact FR II radio source (Brunetti et al. 2001 ) indicates that 
the magnetic field in the radio lobes is very near equipartition; 
equipartitio n cond itions in GPS are also preferred by Sn ellen 
& Schilizzi ( 2002 ). Standard formulae (Pacholczyk 1970 ) were 
used under the assumption that relativistic particles and mag- 
netic fields occupy the same volume (0 = 1), and that the 
amount of energy in heavy particles is the same as that in elec- 
trons (k = 1). We assumed cutoffs in the relativistic particle 
energies at 10 MHz and 100 GHz. The major uncertainty in the 
estimation of the equipartition field strength comes from the 
assumption of the source depth which we assumed to be equal 
to the smaller transverse dimension of the lobe regions where 
we estimated the magnetic field. 

With these assumptions the derived equipartition magnetic 
field is 4.1 milligauss in the east hot spot region, decreasing 
to 3.3 milligauss in the direction of the core. In the less bright 
west hot spot we find 3.5 milligauss and even smaller values 
(2.2-2.7 milligauss) in the extended W lobe regions. The un- 
certainties are mainly due to the non-uniform spectrum present 
in this region. With these results the estimated radiative age of 
older emitting electrons is 3.0 x 10 3 yrs in the east lobe (near 
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Fig. 8. The spectrum of lobes, hot spots and core. The steepen- 
ing of both lobes is clearly visible above 5 GHz. 



the core) and 4.5-4.9 x 10 3 yrs in the northern region of the 
west lobe. An estimate of the uncertainty of these values is not 
obvious since errors are mainly model dependent. 

3.4. Polarization 

After removing the instrumental polarization we did not detect 
any polarized signal from the source. The 3<x upper limit on 
the polarized flux density is < 0.140 mJy/beam at 5 GHz. This 
limit corresponds to < 0.3% polarization at the brighter, eastern 
hot spot, to < 0.6% polarization at the core, and < 3% at typical 
locations in the lobes. 

3.5. Variability 

The sum of model-fit components in the two epochs is 1.21 and 
1 .24 Jy respectively, which gives an insignificant difference of 
just 2%. Given that many of the components are extended, such 
a small variation in the total flux density over 5 years seems 
reasonable. The core component is compact (size ~ 0.81 mas) 
and could be more variable. In fact from the model-fit results 
(Table [l]) the core varies from 20+1 mJy in 1995 to 24 ± 
1 mJy in 2000. Although this variation is only significant at 
the 3<x level, if it is real then it indicates a change by 20%. 
This variation is consistent with the limits on core variability 
implied by the integrated intensity observations of Fassnacht 
& Taylor (2001) from a small sample of CSOs studied over 
the course of eight months. Fassnacht & Taylor found the rms 
variation in the integrated flux density to be 0.7%. Assuming 
all that variation was due to the core component, and given 



Better limits on the variability of core components in CSOs will 
require VLBI monitoring campaigns, rather than low resolution 
studies only sensitive to the integrated flux density. The radio 
galaxy 4C 3 1 .04 is an excellent candidate for such a study since 
the core is comparatively strong and isolated from nearby jet 
components. 

4. Discussion 

4.1. Ages of Compact Symmetric Objects 

We derive the kinematic age of 4C 31.04 to be 550 years. This 
estimate is within the 300 to 1200 year range of ages found in 



the literature (Polatidis & Conway 2002; Owsianik & Conway 



1998; Taylor et al. 2000). It is interesting to ponder the nar- 



rowness of this range. Even younger sources would need to be 
smaller, and/or have faster hot spot advance speeds. Currently 
there is a lower limit on the size of a source of a few mas 
(~ 10 pc) that we can readily classify as a CSO by identifying 
a core component between two hot spots. Assuming that the 
hot spot advance speed through the AGN environment is sim- 
ilar (0.2 c) for an ultra-compact CSO (the denser environment 
on such small scales should make this a conservative assump- 
tion) then we can derive a lower limit to the age of observed 
CSOs of 100 years. There is also an observational upper limit 
on the size of CSOs. This is the size beyond which the CSO 
becomes difficult to observe with VLBI arrays since an inter- 
ferometer is only sensitive to certain spatial frequencies. VLBI 
arrays have difficulty imaging components larger than about 40 
beams, or sources lacking a strong compact component. By the 
time CSOs reach 1 kpc in size the lobes will be largely resolved 
out and the hot spots will be resolved to the extent that a di- 
rect detection is impossible without phase referencing. Velocity 
measurements of heavily resolved components are also much 
more difficult. These parameters give a conservative upperlimit 
on the source age of 10000 years. Future advances in Space 
VLBI may discover radio sources younger than 100 years, and 
the E-MERLIN and EVLA projects may provide sufficient res- 
olution and sensitivity to extend the upper boundary of sources 
for which we can determine a kinematic age. We also note that 
the advance speed of both east and west hot spots in 4C 3 1 .04 
at 0.33 kr} c are at the high end of those measured in CSOs 
to date, corresponding to a hot spot separation velocity of 0.39 
h^L c. This should be confirmed with future epochs. 

It is of interest to compare the dynamic age of 550 years 
- derived from the study of hot spot velocities, to the radiative 
age - which we estimate from synchrotron losses in 4C 31.04 
to be ~ 4000 years. The radiative age is nearly an order of 
magnitude larger than the dynamic age estimate, indicating that 
some of our assumptions are in error. 

First, we note that the kinematic age estimate rests on the 
assumption that the velocity has been constant over the lifetime 
of the source. Although Readhead et al (1996b) argue that only 
mild evolution in velocity is likely, this assumption may not 
hold in 4C 3 1 .04. As discussed in Sect. \.2 , there is evidence 



an average core fraction of 0.03 (Taylor et al 1996), leads to 



an estimate for the core variability in CSOs of less than 23%. 



from the radio morphology of a strong interaction between the 
jet and the ISM. This could be making the jet wander around 
within the lobe, forming new hot spots which drill quickly out 
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for some brief time before being shut off. Further suggestion 
that the kinematic age may be an underestimate comes from the 
west lobe dimension: the weakest detectable emission north- 
ward is about 60 mas which would require a lobe expanding at 
0.41c. 

On the other hand, it is possible that the radiative age is 
overestimated by the equipartition field strength or the value of 
the critical frequency. While the 22 GHz data suggest that the 
value of the critical frequency is quite accurate, we note that 
the equipartition field could easily be higher if either k > 1 
or <p < 1 . In this case the radiative age would be much more 
similar to the kinematic estimate; for instance H eq ~ 6 mil- 
ligauss (given by k — 20 and (f> = 0.7) would yield ~ 1000 
years. Furthermore, this result was obtained assuming no en- 
ergy losses due to Inverse Compton interactions. Even if scat- 
tering with the microwave background can be neglected with 
no harm, it is possible that in such a small source some inter- 
action is present between relativistic electrons in the lobes and 
photons from the active core. The radiative age would then be 
smaller than calculated. 

It is not easy to infer from the present data which estimate 
is the most reliable. We hope to address the question of the dy- 
namic age with additional observations in the next years. Other 
studies of magnetic fields in CSOs would help to understand 
how correct are the assumptions going into the radiative age es- 
timate. Also, high resolution low frequency data could help to 
estimate H from self-absorption analysis, while high frequency 
data will allow a better estimate of the critical frequency. 



4.2. Morphology and interaction with ISM 

The radio morphology of the two lobes is quite different. It 
is likely that their asymmetry is related to a slightly different 
ambient medium. In particular, the interaction between the rel- 
ativistic jet and the ISM can lead to the formation of lobes of 
various morphologies. 

Ac cordi ng to computer simulations presented by Bicknell 
et al. (2002), the extended, low brightness morphology of the 
western lobe could be produced by a jet interacting with many 
ISM clouds. This interaction is not so strong as to destroy the 
jet but would cause: (1) the presence of a hot spot fainter than 
usual and (2) the spread of a large fraction of relativistic plasma 
carried by the jet and therefore the formation of a diffuse, faint 
lobe. 



In the same scenario, the more compact and regular mor- 
phology of the eastern lobe suggests that it is undergoing less 
interaction with the ISM. The deep hole close to the hot spot 
becomes an even more interesting feature. Could it be due to 
foreground absorption or is it truly a large region in the lobe 
devoid of any radio emitting plasma? The steep spectral index 
of 1.7 + 0.04 found from the spectral index image (Fig. ^), ar- 
gues against foreground absorption by thermal electrons which 
would produce an inverted spectrum. This leaves us with a re- 
markable hole in the lobe. It would be interesting to search for 
the presence of dense molecular gas that might be impenetrable 
by the radio plasma. 



4.3. Evolution of CSOs 

Despite uncertainties and some discrepancy between kinematic 
and spectral age for 4C 3 1 .04, our results clearly indicate that 
it is indeed a very young source. Most of the known Compact 
Symmetric Objects have been selected from flux limited sam- 
ples, so that they are usually quite bright and they will likely 
evolve into FRII radio galaxies. It is interesting to ask what 
will be the fate of the weaker 4C 3 1 .04, and try to learn more on 
the relation between 'child' and 'adult' stages of radio sources. 

One of the problems with the youth scenario is given by the 
relative numbers of small sources (CSOs, CSS) and large clas- 
sical double radio galaxies. This has been explained by mod- 



els (e.g. Fanti et al. 1995) in which the total power P of the 



source decreases as a power law function of linear size LS : 
P oc (LSy h . As a consequence of variations of external den- 
sity and pressure, the index h is not very strictly constrained 
and it may actually change on different scales and for different 
sources. 

If we roughly use P oc (LSy 05 as a law for radio power 
evolution (O'Dea |1998| ), we find that 4C 31.04 will drop be- 
low the power of lO 24 5 W/Hz as soon as it will be ~40 times 
larger than its actual size. In other words, by the time 4C 3 1 .04 
will be 5 kpc large, its power will be similar to that of edge 
darkened FRI radio galaxies. Even if this threshold is actually 
not so sharp and the numbers vary, it is interesting to specu- 
late that 4C 31.04 may have an FRI morphology in the future. 
In fact, based on pc-scale jet velocities and properties, G01 ar- 



gued that there are no significant differences between FR I and 
FRII radiogalaxies on pc-scales. It is then interesting to note 
that a possible future FR I radio galaxy now has the same edge 
brightened look of CSOs evolving into FRIIs. Thus, even in a 
completely different evolutionary stage, there would seem to be 
no difference on the small scale between FRIs and FRIIs. The 
discovery of other CSOs among samples of low -power sources 
is needed to better understand the evolution of low -power radio 
galaxies. 

4.4. Lack of Polarization 



Synchrotron emission is intrinsically polarized up to 70%, de- 
pendi ng on the amount of order in the source magnetic fields 
(Burn 1966 ). Why is there so little polariz ed em ission from 
4C 31.04 and other CSOs (Peck & Taylor |2000| )? From the 
Hi observations we know that 4C 31.04 is viewed through a 
moderate density of atomic gas. If there is ionized gas present 
as well along with a magnetic field, then the radiation from 
4C 3 1 .04 that we observe could have been depolarized by a 
change in polarization angle within each 8 MHz IF, or by a 
gradient in the foreground Faraday screen. For bandwidth de- 
polarization to be effective at 5 GHz requires RMs of 10 5 rad 
m~ 2 , while beam-width depolarization requires gradients in the 
RM of 300 rad rrT 2 mas -1 or more. 



5. Conclusions 

From the comparison of the two VLBI images at 5 GHz, a small 
but significant expansion appears to be present in the Compact 
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Symmetric Object 4C 31.04. From the derived expansion rate 
and the current size of the source we calculate an age of 550 
years. A somehow larger value of the age (3000 - 5000 yrs) is 
derived from spectral index analysis. Despite a significant dif- 
ference, both results agree that the source is very young, as al- 
ready demonstrated for other CSOs in the literature. However, 
the relatively high hot spots velocity and its low power make of 
4C 31.04 a peculiar case among CSOs. 

The discrepancy and uncertainty in the age estimates, the 
uncommon presence of a deep 'hole' in the East lobe, the de- 
tection and the variability of the compact core suggest that it 
may be worth observing this source again in the future. New 
VLBI observations will be of interest to look for free-free ab- 
sorption at low frequency, to confirm the expansion, and to look 
for spectral aging at high frequency. 
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